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We present a microscopic treatment of current-induced torques and thermal fluctuations in itin- 
erant ferromagnets based on a functional formulation of the Keldysh formalism. We find that the 
nonequilibrium magnetization dynamics is governed by a stochastic Landau-Lifschitz-Gilbert equa- 
tion with spin transfer torques. We calculate the Gilbert damping parameter a and the non-adiabatic 
spin transfer torque parameter j3 for a model ferromagnet. We find that j3 7^ a, in agreement with 
the results obtained using imaginary-time methods of Kohno, Tatara and Shibata [J. Phys. Soc. 
Japan 75, 113706 (2006)]. We comment on the relationship between s — d and isotropic-Stoner toy 
models of ferromagnet ism and more realistic density-functional-theory models, and on the impli- 
cations of these relationships for predictions of the j3/a ratio which plays a central role in domain 
wall motion. Only for a single-parabolic-band isotropic-Stoner model with an exchange splitting 
that is small compared to the Fermi energy does j3/a approach one. In addition, our microscopic 
formalism incorporates naturally the fluctuations needed in a nonzero-temperature description of 
the magnetization. We find that to first order in the applied electric field, the usual form of thermal 
fluctuations via a phenomenological stochastic magnetic field holds. 

PACS numbers: 72.25.Pn, 72.15.Gd 
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I. INTRODUCTION 

Phenomena related to order-parameter manipulation 
by transport currents have recently received a great 
deal of attention in magnetic metals and magnetic 
semiconductors. Spin transfer torques, which lead to 
current-driven nanomagnet reversal and to domain wall 
motion in narrow wires, have been at the center of this 
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In spin-transfer-torque theory, the usual Landau- 
Lifschitz-Gilbert (LLG) equation of motion for the 
magnetization direction f2 acquires terms corresponding 
to so-called adiabatic and non-adiabatic spin transfer 
torques which are both proportional to current. Both 
torques can be constructed from symmetry arguments by 
requiring that they be orthogonal to the magnetization 
direction and by realizing that the current essentially 
breaks inversion symmetry. The latter implies that, in 
the long- wavelength limit, terms proportional to Vf2 
are allowed in the LLG equation of motion. The adia- 
batic spin transfer torque 3 -^ is defined as — (v s ■ V) Cl, 
where v s is a velocity, proportional to the current, 
that characterizes the efficiency of spin transfer and is 
required for dimensional reasons. The non-adiabatic 
spin transfer torque 1 ^ is given by — /3f2 x (v s • V) Cl and 
is characterized by the dimensionless parameter (3. 

The LLG equation that incorporates both spin transfer 
torques is then given by 



in the long- wavelength low- frequency limit,— where H is 
the effective field and a the Gilbert damping constant. 
The coefficients a and (3 are dissipative, in the sense that 
they relate quantities that are even under time reversal 
to quantities that are odd under time reversal ] 30 : 31 The 
observation that both a and (5 are related to dissipa- 
tion is important in the context of nonzero-temperature 
effects which play an especially important role 2 ^ in ex- 
26,j?c!£ments with magnetic semiconductors ! 27 ' 28 Nonzero 
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temperature is usually accounted for by adding a Gaus- 
sian stochastic magnetic field h to the effective field in 
the LLG equation, 3 ^ 3 ^! 3 ^ 3 !^^ T he strength of this 
noisy magnetic field is related to the Gilbert damping by 
the fluctuation-dissipation theorem, which ensures that 
the system is characterized in equilibrium by the appro- 
priate Boltzmann distribution. A priori the generaliza- 
tion of Eq. ((T|) to nonzero temperatures is not clear, since 
both a and (3 correspond to dissipative processes and 
the fluctuation-dissipation theorem need not be valid be- 
cause the nonzero current implies that the system is out 
of equilibrium. 

As noted in the literature on current-driven domain 
wall motion ; 12 : 17 the case of (3 = a is special because both 
sides of Eq. JT]) then contain the "co-moving" derivative 
D/Dt = d/dt + v s -V, so that the equation of motion ad- 
mits solutions £ld(t) = f2o( x ~ v s0 where fio( x ) is a time- 
independent solution of the LLG equation in the absence 
of currents. The solution Cl^(t) corresponds to "drift" of 
static magnetization textures with velocity v s . Arguing 
that these solutions must exist Barnes and Maekawa 1 ^ 
claim that f3 — a. However, in realistic systems there 
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is no Galilean invariance 4 - that requires the existence of 
such solutions and therefore in general (3 ^ a. 18 ' 19 In- 
stead, the non-adiabatic spin transfer torque acquires 
contributions from all microscopic processes that violate 
spin conservation and therefore correspond to terms in 
the microscopic Hamiltonian that are not invariant un- 
der spin rotations. Such processes also contribute to 
the Gilbert damping term and therefore in principle any 
nonzero Gilbert damping parameter a implies nonzero 
non-adiabatic spin transfer torques, as we show in our 
specific microscopic model calculations. 

The LLG equation [Eq. |T])] is motivated mainly 
by symmetry considerations and contains four different 
quantities whose meaning can be specified precisely only 
by a microscopic theory which details, at least in princi- 
ple, precisely how they should be evaluated given the full 
system Hamiltonian. These quantities are: i) the effec- 
tive magnetic field H, ii) the transport spin velocity v s , 
Hi) the Gilbert damping parameter a, and iv) the non- 
adiabatic spin-transfer torque parameter /?. The effective 
magnetic field H includes the external magnetic field and 
additional contributions due to magnctostatic interac- 
tions and magnetocrystalline anisotropy. The physics of 
H is well understood 4 ^ and not the subject of this paper. 
The three remaining quantities emerge in a microscopic 
theory from the slow (up to first order in time deriva- 
tives or frequency to) smooth (up to first order in space 
derivatives or wavevector q) response of the magnetiza- 
tion direction to an external magnetic field, in the pres- 
ence of an external electric field which drives a transport 
current. The coefficient a then emerges as the ratio of 
the reactive and dissipative contributions that appear at 
first order in lu in this response function. When spin-orbit 
interactions are neglected, it is easy to verify that the co- 
efficient of the reactive term in the total spin-response is 
the unperturbed spin-density, explaining the unit value 
of this coefficient in Eq. |T]). The two first order space 
derivative terms in this equation reflect respectively the 
change in the reactive and dissipative response due to an 
external electric field. Like its zero-current counterpart, 
the current-related reactive terms can be understood in 
quite general terms based only on spin-conservation con- 
siderations, while the dissipative term is sensitive to mi- 
croscopic details. 

As explained above the condition (3 = a corresponds 
to Galilean invariance at a macroscopic level. Since the 
dissipative terms emerge from spin-dependent disorder 
(or spin-independent disorder when spin-orbit interac- 
tions are included in the cyrstal band structure), it is 
clear that Galilean invariance does not hold microscop- 
ically. Our calculations show that (3 « a can occur in 
models with very specific properties 1 ^, but does not oc- 
cur in general. For example (3 ks a occurs in the specific 
toy model that we study below only when the ferromag- 
netism is weak in the sense that the exchange splitting is 
much smaller than the Fermi energyii. We believe that 
we obtain this result only because the model has isotropic 
parabolic bands and that (3 = a (macroscopic Galilean 



invariance) can occur only accidentally in systems with 
either realistic bands or realistic disorder. In the impor- 
tant transition metal ferromagnet spintronic materials in 
particular, we will argue that orbitals which have dom- 
inant d-character contribute more strongly to the mag- 
netization than to transport and that (3 will tend to be 
larger than a as a consequence. 

In this paper we present a microscopic derivation of the 
equation of motion in Eq. |T]) of the direction of magne- 
tization in the presence of current and at finite temper- 
ature. We use a functional formulation of the Keldysh 
non-equilibrium formalis m 41 ! 42 which leads in a natural 
way to the path-integral formulation of stochastic differ- 
ential equations i 43 ' 44 Within our microscopic treatment 
the dissipative nature of a and (3 is explicit because, as 
briefly discussed above and to be shown in more detail, 
they follow from the dissipative part of the spin-density 
spin-density response function and photon two-magnon 
interaction vertex, respectively. We focus on the simple 
microscopic toy model used in previous wor k 17 i 18 that is 
intended to provide a qualitative description of a generic 
ferromagnet and includes disorder and short-range re- 
pulsive electron-electron interactions. The model's fer- 
romagnctism is treated at the level of Stoner mean-field 
theory. For the disorder we use the same model as in 
Ref. [18| and, where applicable, our results for a and 
(3 agree with theirs. Our random phase approximation 
treatment of Stoner quasiparticle fluctuations evinces the 
equivalence of Stoner and s — d models, in the sense that 
in both models the quantities a and (3 are determined by 
the same response function. In particular a ^ (3 for both 
models in general. 

One benefit of the concinnity of the functional formu- 
lation of the Keldysh formalism is that it enables a nat- 
ural determination of the thermal fluctuations without 
explicitly appealing to the fluctuation-dissipation theo- 
rem. We find that to lowest order in the applied electric 
field the form usually assumed for the strength of the 
fluctuations holds and that there is no contribution to 
the white-noise thermal fluctuations that is related to the 
non-adiabatic torque. We emphasize that this formalism, 
which has been reviewed in other publications and ap- 
plied to other problems ; 41 ' 42 is similar in structure to the 
functional formulation of standard equilibrium Green's 
functions for linear response theory, but is more power- 
ful for non-equilibrium and non-linear problems. 

Since the formalism we use may not be familiar to most 
readers, we first present the model and main results in a 
separate section, namely Sec. [TTJ In Sec. IIIII we present 
the formalism and outline the calculations. In the ap- 
pendix we carry out a typical calculation in more detail. 
Both Sec. IIIII and the appendix may be skipped by read- 
ers who are familiar with the formalism or who may be 
more interested in the results obtained. We end in Sec. II VI 
with our conclusions. 
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II. MODEL AND SUMMARY OF RESULTS 



characterized by 



We model the disordered itinerant ferromagnet as elec- 
trons with delta-function-likc repulsive interactions, us- 
ing the Hamiltonian, 



ft 2 V 2 
2m 
1 



-V a (x)T a 



A,, 



-T 2 + V (x) 



t/>(x, t) + -J(x,i) • A(x, t) 



+^(x,t)^](x,t)^(x,t)^ T (x,t)| , (2) 

where for notational convenience we have introduced the 
spinor 



(3) 



In these expressions the Heisenberg operators ip a (x, t) 
annihilate an electron in the spin state labelled by a G 
and obey the usual equal-time commutation re- 
lations. These spin states have their quantization axis 
parallel to an external Zeeman magnetic field in the z- 
direction which contributes A oxt to the energy differ- 
ence between minority and majority spins. Note that in 
Eq. © the Pauli matrices are indicated by r a , and that 
a sum over the repeated index a 6 {x,y,z} is implied. 
The free-electron dispersion at momentum Kk, given by 
£k = h 2 h 2 /2m, is parabolic with an effective mass m (h 
denotes Planck's constant). 

We choose a delta-function interaction with strength 
U because then the field-theoretic procedure to introduce 
the magnetization direction as a dynamic variable is eas- 
ier to implement. This so-called Hubbard-Stratonovich 
transformation^ can also be generalized to spatially non- 
local interactions^ This procedure, to be discussed in 
more detail in the next section, also yields the mean- 
field, i.e., Stoner, saddle-point equation for the exchange- 
interaction contribution to the spin splitting 



A = U 



N F 
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where N-p(x) = [exp (x/k^T) + 1] is the Fermi dis- 
tribution function with k^T the thermal energy, [i is 
the chemical potential that includes a Hartree mean-field 
shift, and p s is the magnetization density. In practice we 
do not explicitly determine the exchange splitting from 
this equation, but simply assume A is a solution whose 
value may be determined from experiment if needed. 
This is another reason for simply using a delta function 
interaction. 

For the disorder we use the same model as Kohno et 
alJ£ in which the spin-dependent disorder potentials are 



K(x)V 6 (x') = <7 <5(x-x')<U 



(5) 



where ■ ■ • indicates averaging over different realizations 
of the disorder. For randomly distributed scatterers 
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a z = n R u: 
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where U{ (u s ) and n\ (n s ) are the strength and density of 
the scatterer charge (spin) component, respectively, and 
S 2 denotes the average scatterer field orientation. Within 
the self-consistent Born approximation the decay rate "f a 
and lifetime r^ c of a plane wave with spin state \a) are 
determined from 
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where the density of states per spin at the Fermi level 
v a = mkF <T /2ir 2 H 2 , and the Fermi wave number kp a — 
v /2m(e F + oM)/h 2 , where M = (A + A cxt )/2 is the total 
spin splitting. 

Finally, the current in our theory is induced by an 
external homogeneous electric field E that, in the London 
gauge, is related to the vector potential by 



A(t) 



cE 
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where u> p is the frequency of the electric field, to be taken 
to zero eventually, and c is the speed of light. In the 
Hamiltonian [Eq. @] the vector potential is minimally 
coupled to the electrons via the charge current-density 
operator 



J(x,t) 
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^t(x,i) V ^(x,t) - (Vv>t( x ,t)) t/>(x,i)]9) 



with — |e| the electron charge. In the above expression 
we have omitted the diamagnctic contribution as it plays 
no role in the following. 

In the next section we derive, starting from the hamil- 
tonian in Eq. ([2]), the equations of motion for long- 
wavelength deviations 5tt of the magnetization direc- 
tion from the collinear ground state, defined by £1 — 
I + 5fl. We find that these transverse deviations obey 
the stochastic equations of motion 
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where a sum over repeated transverse indices a, b £ {x, y} 
is implied and e a f, is the two-dimensional Lcvi-Civita ten- 
sor. The Gilbert damping parameter is given by 
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with the magnetization density p a — A/U. (For the s — 
d model p s corresponds to the carrier spin polarization 
density.) The velocity v s is related to the electric field 
by 



JejE 
mp s 



l T r T 



(12) 



in terms of the density of majority and minority elec- 
trons, denoted by n-f and n±. Using the fact that to 
linear order in the electric field the current densities of 
the majority and minority electrons arc determined from 
j fj = n a \e\ 2 T* c E I 'to, we observe that the expression for v s 
reduces to the usual expression v s = (j-p — j^) / ( — |e| /0 S ) . 
Our result for the /?-parameter reads 







(13) 

Notice that, as expected, only spin-dependent scattering 
contributes to the non-adiabatic torque parameter (3 and 
the Gilbert damping parameter a. 

In addition, we find that the thermal fluctuations via 
the stochastic magnetic field h are determined by 
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(14) 

where the average is over different realizations of the 
noise. We stress that this form for the strength of the 
fluctuations is derived explicitly, without appealing to 
the fluctuation-dissipation theorem. This is important 
since it is not a priori obvious that in the current-carrying 
situation the fluctuation-dissipation theorem holds. The 
form of the strength of the fluctuations in Eq. (TTJJ is 
however of the usual form, i.e., it is of the same form as 
inferred by the equilibrium fluctuation-dissipation theo- 
rem. This result comes about because shot noised con- 
tributions to the magnetization nois o 48 i 49 enter as higher- 
order terms in the applied electric field than the linear 
response in electric field considered here. 

The linear response result in Eq. ()10|) is consistent 
up to 0(6fl) with the Landau-Lifschitz-Gilbert equa- 
tion that includes both non-adiabatic and adiabatic spin 
transfer torques and thermal fluctuations 



( £ + -v s • V 

ot a 



(15) 



where h is a stochastic magnetic field that obeys the 
correlations given by Eq. f| 14[) . 

We end this section by sketching how the various re- 
sults come about. In the theory to be discussed in more 
detail in the next section, the two quantities of inter- 
est are the transverse spin-density spin-density response 
function (or magnon propagator), which determines the 
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FIG. 1: a) Feynman diagram for the transverse spin-density 
spin-density response function. This diagram may equiva- 
lently be thought of as the spin-wave, or magnon, propagator. 
The magnon frequency is denoted by u and its momentum by 
q. b) Photon two-magnon interaction vertex that ultimately 
gives rises to spin transfer torques. (Note that the photons 
correspond to the external electric field used in our theory to 
induce transport current and hence to terms proportional to 
the current within linear response.) This vertex describes the 
interaction of spin waves with frequency uj and momentum 
q with the electric current that is generated by an external 
electric (photon) field of frequency uj p . (The theory of spin 
transfer torques is related to the u> p — > limit of this diagram 
and this limit is taken in the Feynman diagram.) 



Gilbert damping parameter, and the photon two-magnon 
interaction vertex which gives rise to both the adiabatic 
and non-adiabatic spin transfer torques. (Note that the 
photons simply correspond to the external electric field 
in this case.) Feynman diagrams for both of these func- 
tions are given in Fig. [1] (a) and (b). Quite generally, 
the response function in Fig. ([T]) has reactive and dis- 
sipative parts. In the long time and length scale ex- 
pansion corresponding to the LLG equations, the small- 
frequency zero-momentum part of the reactive contribu- 
tion gives rise to the time-derivative on the left-hand side 
of Eq. (TH)]) . (Note that after Fourier transformation fre- 
quencies turn into time derivatives.) The small- frequency 
zero-momentum part of the dissipative contribution to 
the same response function determines the Gilbert damp- 
ing term on the right-hand side of Eq. (fit))) . Physically, 
this dissipative contribution comes from spin waves that 
decay into particle-hole excitations. Energy conservation 
then leads a delta- function- like, i.e., dissipative, contri- 
bution of the form S(huj— £1+62) where fku is the energy of 
the spin wave and t\ — £2 the energy of the particle-hole 
pair. Summing over all possible particle-hole pair en- 
ergies and performing a zero-momentum low-frequency 
expansion then leads to the Gilbert damping term on 
the right-hand side of Eq. (|10j) . Similarly, the spatial 
derivatives on the left-hand side of Eq. (jTUJ) are the re- 
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suit of the reactive contribution to the zero-frequency 
small-momentum behavior of the photon spin-wave in- 
teraction vertex and give rise to the adiabatic spin trans- 
fer torque. The non-adiabatic torque, proportional to (3 
on the right-hand side of Eq. (fT0|) then emerges from the 
dissipative part of the interaction vertex, and gets con- 
tributions from physical processes in which a spin-wave 
interacts with the current and subsequently decays into 
an incoherent particle-hole excitation. 



III. NONEQUILIBRIUM MAGNETIZATION 
DYNAMICS 

In this section we derive the stochastic equation of mo- 
tion for the transverse magnetization in the presence of 
current. We start out by deriving the general equations, 
and subsequently give the results for the long-wavelength 
low- frequency limit. We discuss the equilibrium situa- 
tion, i.e., the case without electric field, and the nonequi- 
librium situation separately. 



A. Stochastic Equations of Motion 

Our starting point is the path-integral expression for 
the coherent-state probability distribution, written as a 
functional integral^ 



P[0*,<M = J ^ T 1^ T ]^r]^l]exp|-5[^,^] 

(16) 

Roughly speaking, this distribution specifies the proba- 
bility for the system to be in the Grassman coherent state 
4>(x, t). The action is expressed in terms of the fermionic 
fields if) and tp* by 



S[ip*,ip] = J dt' Jdx^i/>*{x,t')ih 
-H[i/>*(x,0,</>(x,0] 



, dlp<r(x, t') 



or 



(17) 



The funtional integration in Eq. (fT6|) is over all fields 
evolving forward in time from — oo to t, and back, thereby 
defining the time integration in the action in Eq. (|17|) to 
be over the Keldysh contour C*. 

We rewrite the interaction term as£°- 



(18) 



with h(x, t) an arbitrary unit vector that determines 
the spin quantization axis. Functional integration over 
the latter enforces rotation invariance 5 ^. The interac- 
tion terms on the right-hand side of Eq. (fT5|) are de- 
coupled by writing them as a Gaussian functional inte- 
gral over a density field (p(x, t)) = and spin- 
density field (A(x, t)n(x, if) = U(^*Tip)/2, respec- 
tively. [The precise meaning of the brackets (• ■ •) is 



defined below Eq. (|23|) .] This Hubbard-Stratonovich 
transformatio n 41 ' 45 ! 50 then introduces the density and 
spin density as dynamical variables in the path-integral 
in Eq. (|16p . Density and spin-density amplitude fluctu- 
ations are gapped and can be approximated at low tem- 
peratures and energies by their saddle-point values. For 
the density we then find a Hartrcc-Fock equation, giving 
rise to a mean-field Hartrcc shift which we absorb in the 
chemical potential. For the spin-density amplitude we 
find the saddle-point equation for A in Eq. 

After these steps we ultimately find that the probabil- 
ity distribution is given by 



xcxp<j -S'[^*,i>,h] } , (19) 



where the unit vector Q enters as the boundary condition 
at t' = t on the functional integration over the fluctuat- 
ing magnetization orientation n. We do not explicitly 
indicate the boundary condition on the fermion fields, 
because, as we shall see, the quantity that enters is the 
fermion Green's function which is determined without 
explicitly referring to the boundary conditions. The ac- 
tion S'[ip*, i/', n] is, using the same notation as for the 
hamiltonian in Eq. explicitly given by 



S'[ip*,^,n} = J dt' jdxLp*{x,t') 
+ ^h(x, t') • T + ^r z - y (x) - V a (x)r a 



t d n 2 v 2 
ih- — i 

dt' 2m 



1 



J(x, t')-A(t'] 



<A(x,i') 
(20) 



At this point we note that, if we would add a separate 
Berry-phase term in this action to enforce the angular- 
momentum like quantization of n, the resulting action 
would be the starting point for treating the s — d model. 

We now do perturbation theory around the collinear 
state by writing 

/ Sn x (x,t) 
h(x,t) ~ <5%(x, t) 

V l-\(5n x (x,t)) 2 -\(5n y {x,t)f 

(21) 

and integrate out the electronic fields using second-order 
perturbation theory in Sn a (x, t) , and first-order pertur- 
bation theory in A(x, t). (Note that to find an equation 
of motion for 5n a that is valid up to first-order the action 
needs to be determined up to quadratic terms in 5n a .) 
To this order the effective action for the magnetization 
is written in the form 

S* cff [<5n] - J dtf J d X ^-j Ps Sn 2 a (x,t') 

+ J dt" J dx Sn a (x,t')H a b(x — x.' ;t' ,t")8nij(x ,t") 
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+ / dt" J dx' Sn a (x,t')K ab (x-x';t',t")6n b (x',t") 

(22) 

where a sum over repeated transverse indices a, b S {x, y} 
is again implied. In Eq. (|22|) . the function H a t(x;t, t') is 
determined by the transverse spin-density spin-density 
response function 

U ab (x-x';t,t') = 

^(^(x,t)T a ^(x,t)^(x',t')T b ^(x',t')) , (23) 

shown diagrammatically in Fig.[T](a). In this expression 
the brackets (• ■ •) = Tr [p(— oo) ■ ■ ■] denote an average 
with respect to the density matrix p(— oo) of a system of 
electrons with the action in Eq. ([20]) . with h(x, t) = z and 
A = 0, that is in equilibrium. The function K ab (x;t,t') 
is given by 



(</>t (x, f )r ^(x, t')^(x', t')ni>(x', t')3(x", t")) 



K ab {x - x'; t, t') = J dt" J dx" 
xA ab (x,x',x";t,t',t")-A(t") 



8h 2 c 



(24) 



where the photon two-magnon vertex function 
A ab (x, x',x";t,t',t") 



(25) 



shown Feynman diagram in Fig. [1] (b). 

We note at this stage, that the above procedure, i.e., 
integrating out the fcrmionic degrees of freedom after ex- 
panding the Hubbard-Stratonovich fields using second- 
order perturbation theory, recovers the random phase 
approximation (RPA). The usual structure of the RPA 
response function contains the Stoner enhancement fac- 
tor 1/(11 — U) where II is the zeroth-ordcr "bubble" 
diagram. This form applies to gapped fields such as 
the density-density and longitudinal spin-density spin- 
density response functions, whose fluctuations we have 
neglected. The transverse spin-density spin-density re- 
sponse function does not have this Stoner enhancement 
factor. 

Next, we split the magnetization into semiclassical and 
fluctuating parts, according to 



5n a (x,t±) = Srt a (x,t) ± 



£a(x,t) 



(26) 



where t + and t_ refer to the forward and backward 
branches of the Keldysh contour, respectively. This 
transformation results in the action 



Seff[«U] 



dt' dx 



<5fi a (x, t')£ a (x, t') 



dt' Idx I dt" I dx' \5tl a (x,t') U^ix-x'it' -t'^+K^ix-x'^' -t") &(x',0 



-(-)/ 



— OO 
t 



dt' Idx I dt" ldx'\U^t') Il^ix-x'^' -t") + K [ a l>(x-x';t' -t") SQ b (x',t") 



(+), 



+ / dt' / dx I dt" I dx' \ 2^a(x, t') [n«(x - x'; t' - t") + Kf b {x -x';t'- t") 

I 



6(x',<'' 



(27) 



where the time integrations are now over the real axis perscripts (+) and (— ), respectively, are related to the 

from minus infinity to t. analytic pieces by 

Before we proceed, we make some general statements 

about dealing with functions on the Keldysh contour^ A [± \t,t') = ±0(±(t - t')) [A>{t,t') - A < (t,t')] . (29) 
A general function A{t, t'), with time arguments on the 

Keldysh contour, can be decomposed into its analytic In addition, the Keldysh part, which, as we shall see, 



pieces by means of 

A(t, t') ee 6(t, t')A> (t, t') + 6(t', t)A< (t, t') 



(28) 



determines the strength of the fluctuations, is defined by 
A K {t,t')=[A > (t,t') + A < (t,t')] . (30) 



with 9(t, t') the Heaviside step function on the Keldysh Note that in the effective action [Eq. ([27]) ] the retarded, 
contour. Generally there can also be a piece A s 5(t,t'), advanced, and Keldysh parts of the various functions de- 



but such a general decomposition is not needed here. Re- 
tarded and advanced functions, distinguished by the su- 



pend only on the difference of time arguments (we have 
implicitly taken the limit w p — ► 0). 
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To derive the equation of motion that is obeyed by 
the magnetization 5Q a we perform another Hubbard- 
Stratonovich transformation, and write the part of the 
action that is quadratic in the fluctuations £ a as a Gaus- 
sian functional integral over an auxiliary field rj a which 
will turn out to correspond, up to prefactors, to the 
stochastic magnetic field h. Explicitly, we then have for 
the probability distribution that 



p[<t>\<t>M 



d[8n]d[g\d[rj\ 



xP[ri]e X p{-S e n[5n^,r l } \ ,(31) 



with the effective action 



S eB [Sfl^,rj\ = J dt' y"dx|-^<Jn a (x,0^(x,f) +7 ?a (x,t)Ux,t)j 



— oo 

t 



dt' I dx I dt" I dx' I SQ a (x, t') n^ } (x - x'; t' - t") + {x - x'; t' - t") &(x', t 



(-)/ 



dt' I dx I dt" I dx' \ £ a (x, t') [n<{' (x - x'; t' - t") + K%> (x - x'; t! - *")J «Mx', f' 

I 



-(+)/ 



(32) 



This action is now linear in the fluctuations and the 
funtional integration over these fluctuations leads to a 



constraint that is precisely the equation of motion for 
the magnetization Sfl. We find that 



J 



n 



(+) 



dt 



If 



(-) 

ba 



iV,-i 



dt 



L afc 



-iV,i 



dt 



K 



(-) 

6a 



iV,-i 



dt 



sn b {x,t) 



^5Q a (x,t)=r ]a (x,t) 



(33) 



with n a (,(q, uj) and K ab (q, u>) denoting the Fourier trans- 
forms of n a b(x — x'; t — t') and K ab (x — x';t — t'). respec- 
tively Note that the procedure used in Eq. (|26|) and (|3Tj) . 
which leads ultimately to the above equation of motion, 



circumvents the usual difficulties of deriving an equation 
of motion with dissipative terms from an action. The 
probability distribution for the noise is given by 



P[rj] = cxp 



dt' / dx / dt" / dx' \ 2?y a (x, t') n* 6 (x - x'; t' - t") + Kf b (x - x'; t' - t") r lb (x' , t") 



so that the correlation function of the stochastic magnetic 
field follows as 

(Va(x,t)rib(x' ,t')) noise = - [Uf b (x-x';t-t') 

i 

+K*(x-x';t-t')} . (35) 

The results in Eqs. ([55)) and (|35[) are the main results of 
this subsection. Clearly, our main tasks are now to deter- 



mine the long-wavelength low-frequency behavior of the 
non-equilibrium spin-density spin-density response func- 
tion and the photon two-magnon vertex function. These 
calculations will be outlined in the next two subsections. 
We start out with the equilibrium situation in which the 
electric field is zero and we only need to consider the 
spin-density spin-density response function. 



a) 




b,t' 



FIG. 2: a) Lowest-order diagram, and b) first vertex correc- 
tion to the spin-density spin-density response function. The 
momentum of the spin wave is denoted by q. Impurity scat- 
tering is indicated by the thin dashed line where the x indi- 
cates the position of the impurity. 



B. Equilibrium Situation 



As noted by Kohno et al^L, because the life time of 
the electrons is almost always extremely small compared 
to the spin splitting in metallic ferromagnets, i.e., h^g- <c 
M, we only need to consider the first vertex correction 
to the spin-density spin-density response function. The 
diagrams that contribute to leading order are given in 
Fig. [2] The corresponding expression reads 



n a6 (q;M') 



iA 2 
~8h 



| / ^3 Tr ^G( P + q;t,t'KG(p;t',t)] 
h E ^'/t^w/t^s/ dt "l * , ''Tr[T a G(p 1 +q;^,^'0r a ,G(p 2 + q;^^OT b G(p 2 ;^^n^G(p 1 ;r,^)] 

h a>&{0,x,y,z} J {2lT > J (27r > Jc ™ Jc °° 



(37) 



where the trace is over spin space and tq denotes the 
2x2 identity matrix. The first term in this equation 
corresponds to the lowest-order diagram in Fig. [2] a) and 
the second term to the diagram with the vertex correction 
in Fig. [2] b). 

The Green's function is defined as 

iG ffff /(x-x';t,0 = <W(x,i)C'(x',t')> 
= 6(t, t')Tr [p(-oc# CT (x, t)i>l (x', If) 

-6(t', i)Tr [p(-oo)#, (x', t')4(x, t) 

so that the Fourier transforms of its analytic pieces read 

-iG < (k,u) = A(k,Lu)N F (huj - /i) ; 
iG > (k,uj) = A(k,u)[l-N F (fiw-fi)} , (38) 

where the spectral function A(k, ui) is defined by 

A(k, uj) = i G (+) (k, u) - G (_) (k, w) . (39) 

Finally, the retarded and advanced Green's functions are 
given by 



Hljj ± — £k + Ma ± ih^a 



(40) 



where ftw^ = hw ± iO as usual . 

With these ingredients, the calculation of the retarded, 
advanced, and Keldysh components of the response func- 
tion in Eq. (|36p . is, in principle, straightforward. Some 
details of these calculations are described in the ap- 
pendix. Here we directly present the results. For the 



retarded and advanced components we find that 

4t ) (q^) = n(±)( q ,c) = 



8M 



± in A huj 



MP 



4±)(q, W ) = -n,(t ) (q^) = 



A 2 Ps . 



16M 2 



iftw , 



(41) 



with the Keldysh parts given by 

Il£ e (q ) w)=n£,(q,w) = 

'sj (v 2 + j) + 2Slv m 



inA 2 k n T 



n s u 2 



2M 2 



Ill(q,a;)=l]^(q,w)=0 



(42) 



In order to obtain the Gilbert damping coefficient it is 
sufficient to perform a zero-momentum small-frequency 
expansion of this response function. The first term that 
enters in a long-wavelength expansion is quadratic and 
determines the spin stiffness that is not of interest to use 
here. In addition, in order to determine the fluctuations 
it turns out to be sufficient to obtain the zero-momentum 
zero- frequency part of the Keldysh response function. In- 
serting these results into the full equations of motion in 
Eqs. (|33|) and (|35|) straightforwardly leads to the results 
in Eqs. (|10M14[) . with v s = 0. (In arriving at these fi- 
nal results we have taken the limit A ext <C A.) In the 
next section we consider the situation with an external 
electric field which leads to a nonzero spin-transfer ve- 
locity v s . We end this subsection by noting that, from 
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FIG. 3: Feynman diagrams that contribute to the spin-wave current interaction that gives rise to spin transfer torques, a) 
Lowest-order diagram, and b) and c) leading-order vertex corrections. 



a phcnomcnological viewpoint, Eqs. (fT0|) and (fT4|) arc 
under debat o 36 ' 38 , even for v s = 0. We hope that the 
microscopic derivation presented here sheds new light 
on this controversy. Finally, we note that the tempo- 
ral delta function in Eq. (Q3J) arises by taking the zero- 



frequency limit of the Keldysh part of the spin-density 
spin-density response function. This implies that the 
stochastic magnetic field in Eq. corresponds to a 
Stratonovich stochastic process, rather an Ito oneS. 



C. With Current 

Our next task is to evaluate the function K a b{~x\ t 1 1') 
that is proportional to the photon two-magnon interac- 



tion vertex and hence, from a microscopic point of view, 
ultimately gives rise to spin transfer torques. The rel- 
evant Feynman diagrams are given in Fig. [3] and corre- 
spond to the expression 



dpi e 



Tr [r a G( Pl + q; t, t') Tb G( Pl ; t' , t")G{ Pl ;t", t)} ( Pl • E) 



1 



a' £ {0,ic,y,z} 



dt' 



dt' 



dp 2 

(27T) 



x Tr [r Q G( Pl +q; t, *">a'G(p 2 +q; t"", t') n G(p 2 ; t\ t'")T a ,G( Pl ;l/" \t")G( Pl ;t" , t) ( Pl • E) 
T a G(p 1 +q;t,rOr a ,G(p 2 + q;t''^iOT 6 G(p 2 ;f , ,i'')r a ,G(p 2 ;i'',r)G(p 1 ;i ,/ ',t)(p 2 -E)]| I , (43) 



where the trace is again over spin space. In this ex- 
pression the first, second and third terms correspond to 
the Feynman diagrams in Fig. [3] a), b) and c), respec- 
tively. Determining the low-frequency long-wavelength 
behavior of the retarded, advanced, and Keldysh com- 
ponents from Eq. (|43j) is straightforward but rather te- 
dious. Typical steps in the calculations are illustrated in 
the appendix for the spin-density spin-density response 
function. Here, we directly present the results. Note 
that to obtain the spin transfer torques, and in particu- 
lar the /3-cocfficient that characterizes the non-adiabatic 
spin transfer torque, it is sufficient to perform a zero- 
frequency long- wavelength expansion. 



The results for the various parts of the function 
K a b(m t,t'), that ultimately determine the adiabatic spin 



transfer torque, are given by 



Kf y {q 3 u) = K^,(q ) u) = 
^±)(q,c) = -A'W(q^) 



iA 2 h\e\ (q-E) 
l6mM 2 



(n T rf - n;rf 



(44) 



We note that these off-diagonal parts correspond to the 
reactive part of the photon two-magnon interaction ver- 
tex. The dissipative part that gets contributions from 
decay processes and determines the non-adiabatic spin 
torque is given by 



iA 2 H 2 \e\(q-E 



±- 



16toM 3 
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~ (jnu 2 - n s u 2 s S^ (n T rf u x - n^f i/ T ) ; 
K*(q,w)=K*(q,w)=0. (45) 

From the above expression we note that, to first order 
in the electric field, the Kcldysh part of the photon two- 
magnon interaction vertex is zero. Ultimately, this im- 
plies that the current does not alter the thermal fluctu- 
ations, at least to first order perturbation theory in the 
electric field. Finally, wc remark that inserting the above 
results from Eqs. (|4"4"|) and (|4"S")) into the general equation 
of motion in Eq. (|33[) leads in a straightforward manner 
to the results for v s and (3 presented in Sec. |TTJ 

IV. CONCLUSIONS 

In conclusion, we have presented a general framework 
for the derivation of the effective equations of motion 
for the magnetization direction of a metallic ferromag- 
net, including nonzero-temperature effects and current. 
An important aspect of our approach is that the func- 
tional Kcldysh methods we employ enable us to incorpo- 
rate thermal fluctuations via stochastic forces in a unify- 
ing manner, without explicitly invoking the fluctuation- 
dissipation theorem. As a specific example, we have car- 
ried out detailed calculations for the model of a disor- 
dered itinerant fcrromagnet used by Kohno et alr^. Our 
results for the Gilbert damping parameter and the (3- 
parameter that characterizes the non-adiabatic torque 
are identical to the results found by these authors us- 
ing imaginary-time methods. We have in addition deter- 
mined the thermal fluctuations and found that, although 
the non-adiabatic torque corresponds to a dissipative pro- 
cess and the current-carrying situation makes applica- 
tion of the fluctuation-dissipation theorem questionable, 
to first order in the electric field the usual fluctuation- 
dissipation relation to the Gilbert damping holds. 

The method presented here is quite general, and in the 
near future we intend to apply it also to other models 
of ferromagnets. Wc briefly comment on the generality 
of our results and what to expect for other models. In 
practice ferromagnctism in metals is usually described in 
terms of some combination of ground state and time- 
dependent spin-density-functional (SDF) theory. The 
structure of the ground state theory is then the same as 
that of the saddle-point mean-field equations that arise 
in our theory, with the spin-dependent interaction in our 
theory replaced by the spin-dependence of the exchange- 
correlation potential in SDF theory and our parabolic 
bands replaced by more complex bands specific to a par- 
ticular material. For transition metals the more realistic 
bands of SDF theory are hybridized s and d bands with 
k dependent spin-splitting which tends to be larger in 
bands with dominant d-character. Transition metal fer- 
romagnets are sometimes described by a crude model in 
which hybridization is not explicitly accounted for and 
the d-orbitals arc assumed to be fully spin-polarized. In 
this s — d model the rf-orbitals do not contribute to the 



density of states at the Fermi level since the majority 
spins are fully occupied and the minority spins are empty. 
It follows that the d-orbitals do not contribute to trans- 
port or to any other property that involves only orbitals 
at the Fermi energy. When the formalism of our paper 
is applied to an s — d model rather than to the single- 
band model we discuss, the d orbitals can contribute to 
properties associated with the reactive pieces of the re- 
sponse functions we evaluate, but not to the properties 
that come from the low-energy limits of the dissipative 
response function pieces. The d-orbitals do contribute to 
the coefficient of oj which translates into the LLG pre- 
cessional dynamics time derivative for example. [See the 
last line of Eq. (|4"T]) . This q = time-derivative can be 
interpreted as capturing the Berry phase associated with 
adiabatic spin-dynamics 42.] It follows that the d and s 
orbital contributions to this coefficient are proportional 
to their respective contributions to the total spin-density. 
The d-orbitals of ans-rf model do not, however, con- 
tribute to the reactive adiabatic spin-torque (v s ) term 
because the c?-bands are either full or empty and there- 
fore do not respond to an electric field. 

The a and (3 dissipative parameters are both defined 
as dimensionless ratios of coefficient contributions from 
dissipative and non-dissipative terms in the equation of 
motion Eq. (|4Tj) . Because a parameterizes the ratio of the 
two time-derivative terms it is indirectly altered by the ri- 
bands. In contrast, (3 parameterizes the ratio of the two 
space-derivative terms neither of which has a <i-orbital 
contribution. This is the reason, as noted in previous 
studieS) 10 i 17 ' 18 why j3 tends to be larger than a in s — d 
models, especially when the d orbitals make the domi- 
nant contribution to the spin-density. As we have men- 
tioned previously, and originally shown by Tserkovnyak 
et. ali 17 , for spin-dependent scattering models with 
parabolic dispersion a w (3 in a Stoncr band model when 
the exchange splitting is much smaller than the Fermi 
energy, but not when the d-orbital Berry phase is added 
to the reactive time derivative of an s — d model. It is 
perhaps expected that a should approximately equal to (3 
in this limit since all the ingredients necessary for macro- 
scopic Galilean invariance seem to be present. When the 
spin-polarization is small there is little to distinguish one 
direction of spin-polarization from another and therefore 
one position in a spin texture from another. When the 
bands are parabolic in addition, an external electric field 
simply accelerates the system's center of mass. Explicit 
calculations for the present toy model demonstrate con- 
clusively that the a = (3 condition which corresponds 
to macroscopic Galilean invariance is not satisfied gen- 
erally. For more general spin-dependent disorder models 
or more realistic exchange splitting values a and (3 are 
never equal. 

These considerations to not apply directly to transition 
metal ferromagnets because of s — d hybridization and be- 
cause of the large d-orbital contribution to the minority- 
spin density of states. It is nevertheless true that the two 
reactive term coefficients can be expressed approximately 
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as the sum of s and d orbital contributions. In the ab- 
sence of spin-orbit coupling, the coefficient of dfl/dt in 
Eq. ([1]) is rigorously equal to one because the Berry phase 
is proportional to the total spin-density, including the s 
and the dominant d contribution. Similarly the reactive 
coefficient of Vf2 can be understood^ in terms of the 
cancellation between convective and processional contri- 
butions to spin-dynamics in the static limit. It follows 
that the (i-orbital weight in this reactive coefficient is not 
zero, as in the s — d model, but still relatively smaller than 
the d contribution to the reactive time-derivative coeffi- 
cient. These considerations suggest that (3/ a will tend to 
be larger than one in most transition metal ferromagnets. 
The main challenges to addressing this issue more quan- 
titatively for a specific material are achieving an under- 
standing of the nature of its spin-independent and spin- 
dependent disorder, accounting for the spin-orbit cou- 
pling present in the bands of the perfect crystal, and 
evaluating the vertex corrections (whose essential role is 
established by these toy model calculations) in systems 
with complex band structures. 

I 
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APPENDIX: CALCULATION OF THE 
RETARDED, IMAGINARY AND KELDYSH 
PARTS OF THE SPIN-DENSITY SPIN-DENSITY 
RESPONSE FUNCTION 



The spin-density spin-density response function in Eq. (|36|) is given by 

IU(q; t, «') = n° b (q; t, t') + I^q; t, f) , (A.l) 

where the first term on the right-hand side is the lowest-order diagram in Fig. [2] (a) and the second term is the vertex 
correction in Fig. [2] (b). In the first part of this appendix we evaluate the lowest-order diagram. In the second part 
the vertex correction is calculated. 



1. No vertex corrections 

Without vertex corrections the response function 

?'A 2 f dr> 

n° h (q; t, f) = — J ^ Tr [r a G(p + q; t, t') n G(p; t' , t)] , (A.2) 

is shown diagrammatically in Fig. [2] (a). The analytic pieces are now easily determined from the above result and 
given by 



n af(q;*.*') = ^ / ^Tr[r a G>(p + q:t,Or fc G<(p;t',t)] ; 



iA 2 



dp 



8h J (2tt) 3 



(A.3) 



Using the results in Eqs. (|37H40[) we have for the retarded and advanced components of the zero- momentum Fourier- 
transformed response function that 



ttO,(±)/ \ 



de f de' 



dp 



8 J (2tt) J (2tt) J (2tt) 3 
Expanding for small energies we find that 



N F (e-fx)-N F (e'-n) 



e — e' 



Tr [ Ta A{p, e)T b A(p, e')] 



(A.4) 



II a f> (q, U ) * JJ[L 1^ JJ^ [A r F ( e _ M ) _ NF{e > „)] _Z_Tr [r a A(p, e)r b A(p, e') 
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A 2 fku f de f de' f dp . , , , ., V , s . 



-F 



iA 2 fku f dp 



32tt J (2tt) 



r Tr[r a A(p,^)r & A(p, M )] 



(A.5) 



From this find TL x 'x =H y 'y (q, u>) , and II X '^ ^(q, ui) = — H y 'x (q, w). Carrying out the remaining integrations 

we have, in the limit ^ a /M — ► 0, that 



A 2 p s iA 2 irhiu J n i u "rn + «s« s pi + J + S 2 ^ 

Hjcx (q, - ttt-t- ± 



8M 



A 2 Ps . 



M 2 



16M 2 

The Keldysh component of the response function is in first instance given by 

l2 r de f de' f dp 



n afc (q^) : 



r <5(7iw - e + e') 



4 7 (2?r) 7 (2tt) 7 (2tt) 3 
X {[1 - iV F (e - n)] N F (e' - fx) + N F (e - /x) [1 - JV F (e' - p)]} Tr [r A(p, e)r fc A(p, e')] 

From this we see that n°^ K (q,w) = n°< K (q,w), and IT°> K (q,w) = n°^ K (q,w) = 0. We find that 



n^f(q,o;) 



iirA 2 k B T | nrfv^vy + n s u 2 



S 2 ± [v 2 + , 



M 2 



(A.6) 



(A.7) 



(A. 



Moreover, we have that 



n°, K (q, w) = ±2i [2JV B (M + 1] Imn^ ± )(q, 



(A.9) 



where N B (x) is the Bose distribution function. This is the fluctuation-dissipation theorem, which emerges naturally 
from the formalism. 



2. Vertex correction 

The first-order vertex correction is shown in Fig. [2] (b), and given by 

,2 



a'£{0,x,y,z\ 



dt" 



xTr [r Q G(q + Pl ; t, t'")r a ,G(q + p 2 ; t'" , t') n G(p 2 ;t' , t")r a ,G( Pl ; t" , t)} . (A. 10) 

Before we proceed, we state some rules for calculus involving the Keldysh contour. From an equation like 



A(t,t')=j dt"B{t,t")C{t",t') 



(All) 



we find the analytic pieces as 

/OO p oo 

dt"B ( - + )(t,t")C^(t",t')+ dt"B^{t,t")C ( -\t",t') , (A.12) 
-oo J — oo 

where the retarded and advanced components are defined in Eq. (|29j) . It is then also straightforward to show that 

pOO 

AW{t,t')= dt"B (± \t,t")C [±) {t",t') . (A.13) 
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Using these rules we find from Eq. (|A.10[) 

n i,(±), ST ( J±- f de ' f dpi f dp2 1 

06 lq ' Wj ~ 8 ^'7 {2 n)J {2k) J (2k) 3 J (2k) 3 fiu,± - c + 

x Tr { r a G (+) (q + pi ; e)r a - G> (q + p 2 ; e) + r a G> (q + pi ; e)r ' G (_) (q + p 2 ; e) 
r b G« (p 2 ; e')T a ,G < ( Pl ; e') + r b G < (p 2 ; e'K'G*->(pi; e') 

- T a GW(q + pi; e)T /G<(q + p 2 ; e) + r Q G<(q + pi;e)T a ,G^(q + p 2 ; e) 
r b G (+ )(p 2 ; e')vG > ( Pl ;6') + r b G > (p 2 ; e')r a ,G^(p 1 ;e')\ } . 
The Kcldysh component reads 

de f de' f dpi f dp- 



, K . , IkA 2 



n^(q^) 



a'e{0,x,y,z} 



: S(hoj -e + e') 



(2k) J (2k) J (2k) 3 J (2k) 3 
xTr { [r.G<+)(q + Pi; e)vG > (q + p 2 ; e) + r Q G > (q + Pl ; e)r a ,G«(q + p 2 ; e) 
T b G(+> (p 2 ; e')r , G< ( Pl ; e') + r b G< (p 2 ; e')r a - GH ( P i ; e') 
- t q G (+) (q + pi ; e)T a , G< (q + p 2 ; e) + r a G< (q + pi ; e)T a 'G (_) (q + p 2 ; e) 
T b G (+) (p 2 ;e')T a -G > (p 1 ; e ')+T b G > (p 2 ; e 'K'G ( - ) (pi;e')]} • 



(A.14) 



(A.15) 



We focus now on the imaginary part of ni'j (q, oS) = Hy'y (q, u>) since these determine the Gilbert damping constant. 
Carrying out the momentum and energy integrals results in 



n 



i,(±) 



(q,w) = =F 



iKA 2 hu 



The corresponding Keldysh part is given by 



n^ K (q,^) 



■iA 2 k B T 



n i u 2 v ] v i - n s ulS 2 z v^vi 
M 2 



M 2 



(A.16) 



(A.17) 



Adding the results from Eqs. (|A.6|) and (|A.16[) we get the result presented in Eq. (|Hj) . Similarly, adding the result in 
Eq. fOll to Eq. fP7|> reproduces Eq. (JUJ). ^ 
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